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Abstract 

Quantum interference effects and decoherence mechanisms in single- molecule junctions are ana- 
lyzed employing a nonequilibrium Green's function approach. Electrons tunneling through quasi- 
degenerate states of a nanoscale molecular junction exhibit interference effects. We show that 
electronic-vibrational coupling, inherent to any molecular junction, strongly quenches such in- 
terference effects. As a result, the electrical current can be significantly larger than without 
electronic-vibrational coupling. The analysis reveals that the quenching of quantum interference is 
particularly pronounced if the junction is vibrationally highly excited, e.g. due to current-induced 
nonequilibrium effects in the resonant transport regime. 

PACS numbers: 73.23.-b,85.65.+h,71.38.-k 
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Quantum interference is at the heart of quantum mechanics. While this classically 
counter- intuitive phenomenon has been verified long ago, e.g. by double-slit experiments 
, quantum interference effects recently attracted much attention in electron transport 
through nanostructures, such as, e.g., quantum dots 3)- 7] or single- molecule junctions E - 
When a molecule is contacted by two electrodes, forming a single-molecule junction 121. Il3|. 
two different physical regimes become interconnected: The microscopic realm of a single 
molecule, governed by coherent quantum dynamics, and macroscopic reservoirs of electrons, 
wherein electronic coherences decay rapidly. Hence, the question arises to which extent a 
single- molecule junction preserves quantum coherence and which decoherence mechanisms 
are active. Understanding these mechanisms is crucial for the design of nanoelectronic de- 
vices. 

This question can be addressed by analyzing quantum interference effects that occur in 
a single- molecule junction. Consider e.g. the model molecular junction depicted in Fig. [Jji. 
This junction comprises two electronic states, which are located at the molecular bridge. 
One of them is gerade with respect to the L -H- R symmetry of the junction, the other 
is ungerade. If these states are quasi-degenerate, they provide, similar to a double-slit 
experiment, two different pathways for an electron tunneling through the junction. As the 
respective wave-functions for the outgoing electron thus differ by sign, they destructively 
interfere with each other. As a result, the corresponding tunnel current is suppressed, 
or may even completely vanish {vide infra). In this letter we show that such quantum 
interference effects are quenched by vibrationally-assisted tunneling processes (cf . Fig. [lb and 
db), during which the tunneling electron excites or deexcites vibrational degrees of freedom 
of the junction. Due to the small size and mass of a molecular conductor, vibrational and 
electronic degrees of freedom are typically strongly coupled, resulting in inelastic tunneling 



processes, which strongly influence the respective transport characteristics |13l420l|. This 
aspect distinguishes nanoscale molecular conductors from mesoscopic systems and quantum 
dots. 

We describe vibrationally coupled electron transport through a single molecule (M) that 
is covalently bound to two metal leads (L,R) by a generic model with Hamiltonian H = 
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FIG. 1: (Color online) a): Scheme of a molecular junction at a finite bias voltage $ that exhibits 
destructive interference. The gray shaded areas depict the continuum of occupied electronic states 
in the left (L) and the right (R) lead. Electronic states located at the molecular bridge (M) are 
represented by black bars. The blue and red arrows depict the coupling of these states to the leads, 
which due to the different symmetry of the corresponding orbitals, gerade and ungerade, differ 
by sign (example orbitals are given in the insets), b) and c): Example processes for vibrational 
excitation (red wiggly line) and deexcitation (blue wiggly line) upon inelastic tunneling of an 
electron through the molecular junction. 
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The electronic part of the Hamiltonian, H e \, includes electronic states with energies 6j located 
at the molecular bridge (cf. Fig. UK). These states are coupled by interaction matrix elements 
Vki to electronic states with energies in the leads. Thereby, the operators cj/cj (cyck) 
represent the respective creation/annihilation operators for the states of the molecular bridge 
(the leads). The vibrational degrees of freedom of the junction are described as harmonic 
oscillators with frequencies Q a , and corresponding creation/annihilation operators a' a /a a . 
Here, \ a denotes the electronic- vibrational (vibronic) coupling strength between vibrational 
mode a and the ith state of the bridge. 



To study quantum interference and decoherence effects in molecu 



ploy a Non-Equilibrium Green's Function approach (NEGF) 14H17I. |21| . NEGF theory 



ar junctions, we em- 



allows to describe quasi-degenerate molecular levels, which is crucial for the analysis of the 



respective quantum interference effects. The method is based on (self-consistent) second- 
order perturbation theory in the molecule-lead coupling, and treats electronic-vibrational 
coupling non-perturbatively. The details of our NEGF approach have been outlined pre- 



viously 
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211 ] . Briefly, the approach is based on the small polaron transformation 



21] . The accordingly transformed Hamiltonian H comprises a polaron-shift of 
the electronic states of the molecular bridge, e« = q — X] Q (^L/^a)> a molecule-lead cou- 
pling term, which is renormalized by the shift operators X t = exp(^ Q ,(Aj a /f2 Q )(a Q , — ajj), 
and Hubbard-like electron-electron interaction terms, (\i a \j a /VL a )c\ciC^Cj, but no direct 
electronic-vibrational coupling terms. The vibrationally induced electron-electron interac- 
tions are treated using a non-perturbative approximate scheme (jj]]. In the (anti-)adiabatic 
regime of a molecular junction, the electronic Green's function matrix G can be described as 
Gy(r,r') « G C!ij (r,r')(T c X i (r)X}(T')}. The matrix G c ^(r,r') = -i(T c c i (<r)ct(T')) is de- 
termined by the self-energy matrix S L/Riij (r, r') = EfceL/R V ki V kj9k(r, r / )(T c X j (r')X/(r)), 
where gk denotes the free Green's function of lead state k and T c is the time-ordering opera- 
tor on the Keldysh contour. The correlation functions of the shift-operators Xi are evaluated 
using a cumulant expansion in the dimensionless coupling parameters Xi a /Q a , which in turn 
requires the electronic Green's functions G c ^j. Therefore, we employ a self-consistent solu- 



tion scheme 
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211 ] . This scheme accounts for the molecule- lead coupling in second-order 



(self-consistent) perturbation theory. It is thus exact for vanishing vibronic coupling, and 
due to the small polaron transformation, also for vanishing molecule-lead coupling. The 
electrical current through the junction is given by / = 2e f (de/(27r)) tr-fE^G^ — ^lG^} 



HQ, 



2l|. 



First, we consider a generic model system similar to the one depicted in Fig. with pa- 
rameters that are typical for single-molecule junctions. The model comprises two electronic 
states located at Z\ = 0.404 eV and e 2 = 0.4 eV above the Fermi-level ep = 0. State 1 of 
our model system is coupled to the vibrational mode with coupling strength Ai = 0.06 eV. 
For state 2, we consider A2 = 0. This choice simplifies the analysis of the corresponding 
decoherence mechanism, which occurs in general for Ai 7^ A 2 . These states are coupled to 
the leads in the same way as shown in Fig. [T^,, i.e. with ^l,i/2 = ^r,i = —^r,i = v = 0.1 eV. 
The leads are modelled as semi-infinite tight-binding chains with level-width functions 
r L /R,ij(e) = (^L/R,i^L/R,j/7 2 )Re [y^ 2 - (e - /x l /r) 2 J and 7 = 3eV. The bias voltage $ 
is assumed to symmetrically drop at the contacts, /il = — /iR = $/2. 
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FIG. 2: (Color online) Current-voltage characteristics for a generic model molecular junction that 
exhibits two destructively interfering electronic states (see Fig. [1]) coupled to a single vibrational 
mode. The inset shows the corresponding vibrational excitation characteristics, {a^a). 

Current- voltage characteristics for this model system are shown in Fig. [2j Thereby, the 
solid gray line shows the electrical current obtained for static nuclei ('electronic current'), i.e. 
for the polaron-shifted electronic states with Ai = 0. The comparatively small value of the 
electronic current (0.26 /xA for $ > 2ei/2) is a result of destructive interference effects in this 
molecular junction, which arise due to the symmetry of the molecule- lead coupling and the 
quasidegeneracy of the two states (ei — e 2 < T). If the partial currents corresponding to the 
different pathways through the molecular junction are added incoherently (dashed gray line), 
a much larger value (3.2 /xA for $ > 2E1/2) is obtained. This demonstrates the presence of 
destructive quantum interference, which reduces the current through this junction by more 
than an order of magnitude. 

If we take into account electronic- vibrational coupling, Ai = 0.06 eV, vibrational exci- 
tation and deexcitation processes (cf. Fig. [TJd and [Tfc) contribute to the electron transport 



through the junction 



15 



15KL7I]. These inelastic processes result in an highly excited 



state of the vibrational mode for $ > 2ei (cf. the vibrational excitation characteristics shown 
in the inset of Fig. EJ. The respective current- voltage characteristics, which is depicted by 
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FIG. 3: (Color online) Transmission functions t(e) (dashed black line) corresponding to the currents 
depicted by the solid gray, blue and red lines in Fig. [2] at a high bias voltage. For clarity, the 
graphs for these three different cases are separated by offsets. In addition, the incoherent sum 
of transmission channels t snm (e) (orange lines) and the interference term tmt(e) (green lines) are 
depicted (cf. Eqs. ([3])). 



the solid black line in Fig. [21 exhibits several steps at $ = 2 



6i -\-nQ), n G N, corresponding 



to vibrationally-assisted electron transport processes |14Hl7j|. Intriguingly, the current that 
is obtained including electronic- vibrational coupling (black line), is much larger than the 
electronic current (gray line) and approaches the value obtained by incoherent summation 
of the different electronic pathways (dashed gray line). Electronic- vibrational coupling thus 
results in a complete quenching of destructive interference effects in this model molecular 
junction, or equivalently, leads to strong decoherence. While the results discussed so far 
have been obtained using the full current-induced nonequilibrium vibrational distribution, 
the basic mechanism of vibronic decoherence can also be described within the often em- 
ployed approximation, where the state of the vibrational mode is restricted to its thermal 
equilibrium state after each electron transmission event. The solid blue and red lines in 
Fig. [2] show results employing this simpler approach obtained for temperatures of 10 K and 
5000 K I22J, respectively. 
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To analyze these effects in more detail we employ a thermal equilibrium state for the 
vibration and consider bias voltages $ 3> ei/2, for which the current can be expressed in 
a form similar to Landauer theory, / 2e J^ o (de/27r) i(e), with the transmission function 
t(e) = ztr{TLG > }. This transmission function t(e) can be split into an incoherent sum of 
transmission channels, t sura (e) = iYli^u^^Hi "where is evaluated disregarding the off- 
diagonal elements of Sl/r, and an interference term, tmt(e) = t(e) — t 8wn (e). The interference 
term explicitly shows the interference effects between individual transmission channels en- 
coded in the off-diagonal elements of Sl/r- F° r our specific model system, in the wide-band 
approximation with T = 2z/ 2 /7, we obtain 

*Bum(e) = (r/|e-e 2 + ir|) 2 (3) 

oo 

+A Ii(x)eP m ' 2 (r/ \e-e 1 -lQ + iT\) 2 , 

l=—oo 

t int (e) = ^A^ReKe-^+iTy^e-^-iT)- 1 ]. 

Thereby, the prefactor A = e~( A i/ n2 ^ 27Vvib+1 ) is determined by the average vibrational ex- 
citation iV v ib = (e^ — l)^ 1 and the inverse temperature /3 = (ksT)" 1 , while h{x) = 
Ii{2(\\/Vt 2 ) N vi -b(N vi -b + 1)) denotes the Ith modified Bessel function of the first kind. 

Fig. [3] shows the transmission functions (dashed black lines) corresponding to the currents 
represented by the gray, the solid blue and red line in Fig. [2j In addition to the transmission 
functions, the incoherent sum of transmission channels (orange line) and the interference 
terms (green line) are depicted. Two important observations can be made. First, for non- 
vanishing vibronic coupling, the transmission peak associated with state 1 is fragmented 
into vibrational side-peaks. These side-peaks, which are associated with vibronic tunneling 
processes (Fig. [lb and [lb), appear in the incoherent sum of transmission channels, but 
not in the interference term. Thus, these side-peaks can directly contribute to the total 
transmission probability. Second, the prefactor A = e~( A i/ Q2 ^ 2Arvib+1 ) reduces the main peak 
at e ~ ei/2 in both parts of the transmission function. While for iV v ib — > oo the elastic 
peak in the incoherent sum of the transmission channels is thus halved, the peak in the 
interference term completely vanishes for a high level of vibrational excitation. As a result, 
quantum interference effects are greatly reduced, if the vibration acquires a highly excited 
state, in particular in the resonant tunneling regime (cf. the inset of Fig. [5]). 

In the model considered so far, the decoherence mechanism is based on the coupling of 
the electronic states of the molecular junction to its internal vibrational degrees of freedom. 
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In many cases 23|, |24|, however, decoherence (or dephasing) is a result of external or en- 



vironmental degrees of freedom. Such a mechanism can be described by a different model, 
where the electronic states are coupled to a bath of harmonic oscillators, mimicking, e.g., the 
coupling to the phonons in the leads [25( or a surrounding solvent in an electrochemically 



gated molecular junction 25|, |26|. The current- volt age characteristics obtained for such a 



model is shown in in Fig. [2] by the dashed green line. Thereby, the overall vibronic coupling 
strength (as specified by the reorganization energy, 0.36 meV) is the same as for the single- 
mode model considered above but has been equally distributed over 100 vibrational modes, 
which frequencies range from 1 meV to 250 meV with an equidistant spacing of 2.5 meV. The 
comparison to the corresponding single-mode case (solid blue line) shows that the coupling 
to a multitude of vibrational modes induces even stronger decoherence effects. Within the 
analysis of interference discussed above, this mechanism is encoded in the prefactor A of Eqs. 
03]), which for the present model reads A = e~^^ x ^^ n ° ! ^ 2Nvib ' a+1 \ Due to the low-frequency 
modes, it is considerably smaller than for the single-mode model. 

The results discussed so far for a generic model system show that electronic-vibrational 
coupling leads to decoherence in single-molecule junctions and elucidates the underlying 
mechanisms. As a specific example for a single-molecule junction that exhibits quantum 
interference and decoherence, we consider a o-biphenylacetylenedithiolate molecule bound 
to gold electrodes (see Fig. H|). The model parameters of this junction were determined by 
first-principles electronic structure calculations, where a detailed description of our method- 

n 

ology can be found in Ref. [27j. The junction exhibits 36 active vibrational modes coupled 
to a multitude of closely-lying electronic states that interfere with each other. Test calcu- 
lations show that for bias voltages within a range of |$| < 3 V the three electronic states 
corresponding to the molecular orbitals depicted in Fig. H^) are sufficient to describe the 
electronic current in this junction. The comparison of the electronic current obtained within 
this three-state model (solid gray line) with the significantly larger current that is obtained 
for the transport through the individual states (dashed lines) demonstrates the presence 
of pronounced interference effects among the three pathways in this junction. Analysis 
shows that the destructive interference of the highest occupied molecular orbital (HOMO) 
and HOMO-1 are due to their quasidegeneracy and the fact that their wavefunctions are 
symmetric and antisymmetric combinations, respectively. In addition, the electronic state 
corresponding to HOMO-4 exhibits destructive interference with the other two orbitals, due 
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FIG. 4: (Color online) Electron Transport in an o-biphenylacetylenedithiolate molecular junction, 
a): relevant molecular orbitals, where the blue and the red color represent different signs, b): 
current-voltage characteristics for different scenarios as described in the legend. 

to the broadening of the molecular levels. If we account for electronic-vibrational coupling, 
the current-voltage characteristics represented by the solid blue and turquoise line in Fig. 
H] are obtained, where all 36 vibrational degrees of freedom are evaluated in thermal equi- 
librium at 10 K and 300 K, respectively. The associated current- voltage characteristics are 
significantly larger than the electronic current. Moreover, for an increasing temperature of 
the vibrational modes, the current through the junction approaches the significantly higher 
values obtained for the individual channels, given by the dashed lines. As in the generic 
model system, vibronic interactions thus result in a strong decoherence of interfering elec- 
tronic states, which is greatly enhanced by the temperature of the modes or equivalently 
their level of vibrational excitation. 

In summary, single- molecule junctions often exhibit multiple closely-lying electronic states 
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that provide interfering pathways for the tunneling of electrons. Our analysis, obtained for 
both a generic model and a specific example, show that electronic-vibrational coupling leads 
to decoherence and quenching of these electronic quantum interference effects. As a re- 
sult, the electrical current in molecular junctions can be significantly larger than without 
electronic-vibrational coupling. This effect is particularly pronounced if the junction is vi- 
brationally highly excited, e.g. due to current-induced nonequilibrium effects in the resonant 
transport regime. It should also be emphasized that this decoherence mechanism is an in- 
trinsic property of a single-molecule junction. 

We gratefully acknowledge fruitful discussions with S. Ballmann, B. Kubala and P. B. 
Coto, and generous allocation of computing time by the Leibniz Rechenzentrum Munchen 
(LRZ). The work at the Friedrich-Alexander Universitat Erlangen-Ntirnberg was carried out 
in the framework of the Cluster of Excellence " Engineering of Advanced Materials" . 
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